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ABSTRACT 

We use the LUQAS sample (Kim et al. 2004), a set of 27 high- resolution 
and high signal-to-noise QSO absorption spectra at a median redshift of 
z = 2.25, and the data from Croft et al. (2002) at a median redshift of 
z = 2.72, together with a large suite of high- resolution large box-size hydro- 
dynamical simulations, to estimate the linear dark matter power spectrum on 
scales 0.003 s/km < k < 0.03 s/km. Our re-analysis of the Croft et al. data 
agrees well with their results if we assume the same mean optical depth and 
gas temperature-density relation. The inferred linear dark matter power spec- 
trum at 2; = 2.72 also agrees with that inferred from LUQAS at lower redshift 
if we assume that the increase of the amplitude is due to gravitational growth 
between these redshifts. We further argue that the smaller mean optical depth 
measured from high-resolution spectra is more accurate than the larger value 
obtained from low-resolution spectra by Press et al. (1993) which Croft et 
al. used. For the smaller optical depth we obtain a « 20% higher value for the 
rms fluctuation amplitude of the matter density. By combining the amplitude 
of the matter power spectrum inferred from the Lya forest with the amplitude 
on large scales inferred from measurements of the CMB we obtain constraints 
on the primordial spectral index n and the normalisation erg. For values of the 
mean optical depth favoured by high-resolution spectra, the inferred linear 
power spectrum is consistent with a ACDM model with a scale-free (n = 1) 
primordial power spectrum. 

Key words: Cosmology: intergalactic medium - large-scale structure of uni- 
verse - quasars: absorption lines 



1 INTRODUCTION 

The prominent absorption features blue-ward of the 
Lyman-a emission in the spectra of high-redshift 
quasars (QSOs) are now generally believed to arise from 
smooth density fluctuations of a photoionised warm in- 
tergalactic medium (see Ranch 1998 and Weinberg et 
al. 1999 for reviews). This has opened up the possibil- 
ity to probe the density fluctuation of matter with the 
flux power spectrum of QSO absorption lines (Croft et 
al. 1998, Hui 1999, Croft et al. 1999b, McDonald et al. 
2000 [MOO], Hui et al. 2001, Croft et al. 2002 [C02], 
McDonald 2003, Viel et al. 2003). 

The flux power spectrum is mainly sensitive to the 



slope and amplitude of the linear dark matter power 
spectrum for wave-numbers in the range 0.002 s/km < 
k < 0.05 s/km. Croft et al. (1999) inferred an amplitude 
and slope which was consistent with a COBE normalised 
ACDM model with a primordial scale invariant fluctu- 
ation spectrum (Phillips et al. 2001). MOO and C02, 
using a larger sample of better quality data, found a 
somewhat shallower slope and smaller fluctuation am- 
plitude. The WMAP team used the later data in com- 
bination with their CMB data to claim that there is 
evidence for a tilted primordial CMB-normalised fluc- 
tuation spectrum {n < 1) and/or a running spectral 
index (Bennet et al. 2003; Spergel et al. 2003; Verde et 
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al. 2003) . A number of authors have argued that the er- 
rors in the inferred dark matter (DM) power spectrum 
have been underestimated (Zaldarriaga, Scoccimarro & 
Hui, 2003; Zaldarriaga, Hui & Tegmark, 2001; Gnedin 
& Hamilton 2002; Seljak, McDonald & Makarov 2003). 
We here use a suite of high-resolution hydro-dynamical 
simulations and the flux power spectrum obtained from 
LUQAS (Large Sample of UVES QSO Absorption Spec- 
tra), together with the published flux power spectrum 
of C02 to further investigate these issues. 

The plan of the paper is as follows. In Section 2 
we describe the two data sets and give an overview of 
the determinations of the effective optical depth found 
in the literature. The flux power spectra obtained from 
the hydro-dynamical simulations are discussed in Sec- 
tion 3. Section 4 describes the method and uncertain- 
ties of inferring the linear matter power spectrum. In 
Section 5 we present our results, compare with previous 
results and discuss implications for as and n. Section 6 
contains a summary and our conclusions. 
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2 THE OBSERVED FLUX POWER 

SPECTRUM AND EFFECTIVE OPTICAL 
DEPTH 

We will use estimates of the flux power spectrum from 
two different data sets: the LUQAS sample of VLT spec- 
tra compiled by K04 and the sample of Keck spectra 
compiled by C02. Below, we will describe both samples 
in more detail. 

2.1 The LUQAS sample 

LUQAS (Large Sample of UVES Quasar Absorption 
Spectra) compiled by K04, consists of 27 high-resolution 
spectra taken with the VLT UVES spectrograph. The 
median redshift of the Lyman-Q forest probed by the 
spectra is {z) — 2.25 and the total redshift path cov- 
ered is Az — 13.75. The S/N varies as a function of 
wavelength, but in the forest region it is usually larger 
than 50. The flux power spectrum of the LUQAS sam- 
ple has been analysed in K04, while a study of the flux 
bispectrum can be found in Viel et al. (2004b). 

The LUQAS sample probes somewhat smaller red- 
shifts than the C02 sample, but, as discussed by K04, 
in the redshift range where the two samples overlap, the 
estimated flux power spectra agree well in the range of 
wavenumbers not affected by differences in resolution 
and S/N. Here, we will use a subset of the LUQAS sam- 
ple for which we have selected all QSO spectra regions 
in the redshift range 2 < z < 2.3, to maximise the 
contrast in redshift to the C02 sample and to further 
investigate the redshift evolution. 16 QSOs contribute 
to this subsample (Figure 1 in K04) and the median 
redshift is (z) = 2.125. In Table 1, we give the 3D flux 
power spectrum Ap{k) = P|°(fc) P/{2tv^) for this sam- 
ple, using the flux estimator F/ (F) — 1 (denoted F2 in 
K04), where F is the flux of the continuum-fitted spec- 
tra. The 3D flux power spectrum is obtained from the 
ID flux power spectrum by differentiation in the usual 
way, viz. 



Figure 1. Triangles show the effective optical depth r^g = 
- In {F) for the LUQAS sample (also given in Table 2). The 
dashed curve is the result of Kim et al. (2002) while the long- 
dashed and dot-dashed curves show the result of Schaye et al. 
(2003), with and without the removal of pixels contaminated 
by metals, respectively. The dotted curve is the effective op- 
tical depth obtained by Press, Rybicki &; Schneider (1993). 
The solid curve shows the results from Bernardi et al. (2003). 
The empty diamond is the result of Tytler et al. (2004). The 
vertical extent of the two shaded regions indicates the values 
used in our analysis. 



Note that peculiar velocities and thermal broad- 
ening make the flux field anisotropic and that Pp^ is 
thus not the true 3D power spectrum of the flux. The 
flux power spectrum has been calculated for the same 
wavenumbers as in C02. 

2.2 The Croft et al. sample 

The Croft et al. (2002) sample (C02) consists of 30 
Keck HIRES spectra and 23 Keck LRIS spectra. We will 
here use the flux power spectrum which C02 derived for 
what they call their 'fiducial' sample. This sample has 
a median redshift (z) = 2.72, spans the redshift range 

2.3 < z < 3.2, and has a total redshift path of Az — 25. 
For more details see C02. 

2.3 Statistical and systematic errors in the 
observed flux power spectrum 

Unless stated otherwise, we estimate our statistical er- 
rors with a jack-knife estimator. The main systematic 
errors affecting estimates of the flux power spectrum 
due to absorption by the Lyman-a forest are contin- 
uum fitting and the presence of metal lines and damped 
Lyman-Qf systems. These effects have been investigated 
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Table 1. The ID and 3D power spectrum of the flux F = 
exp(— r)/ (exp(— t)) — 1 of the LUQAS sample at 2 = 2.125. 
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Table 2. Mean flux and effective optical depth of the LUQAS 
sample. 



redshift range 




in 




2.0 < 2 < 2.3 


2.125 


0.849 ± 0.008 


0.163 ±0.009 


2.3 < z < 2.6 


2.44 


0.800 ± 0.008 


0.223 ±0.014 


2.55 < z < 3.0 


2.72 


0.730 ±0.011 


0.315 ±0.015 



by e.g. C02 and K04. The main conclusions of K04 arc 
as follows: i) continuum fitting uncertainties for high- 
resolution Echelle spectra strongly affect the flux power 
spectrum at scales k < 0.003 s/km; ii) the contribution 
of metal lines is less than 10% at scales k < 0.01 s/km 
but rises significantly (up to 50%) at smaller scales; 
Hi) damped Lyman-a systems appear to only mildly 
affect the estimate of the flux power spectrum. In or- 
der to minimise uncertainties due to continuum fitting 
and metal lines, and to avoid dealing with the problem- 
atic thermal cut-off at small scales, wc will only use the 
range of wavenumbers 0.003 < k (s/km) < 0.03 for our 
analysis. 

2.4 The observed eflfective optical depth 

As pointed out by Croft et al. (1998), C02 and Seljak, 
McDonald & Makarov (2003) and discussed in detail in 
Sections 4 and 4.2, the assumed effective optical depth, 
Teff = — in (F) , has a large influence on the amplitude of 
the inferred matter power spectrum. The observed val- 
ues of TefT have largo statistical and probably also not 
yet fully understood systematic errors. The main uncer- 
tainty in determining the effective optical depth comes 
from the continuum fitting procedure and the Poisson 
noise due to the large variations from line-of-sight to 
hne-of-sight (Zuo & Bond 1994, Viel et al. 2004a, Tytler 
et al. 2004). 

In Figure 1, we show results for the redshift evo- 
lution of Teff from a number of observational studies in 



Table 3. Grid of cosmological simulation parameters. 
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the literature. The dotted curve is from a large sample 
of low-resolution spectra compiled by Press, Rybicki & 
Schneider (1993, PRS), while the solid line is from a 
sample of low resolution spectra drawn from SDSS by 
Bernardi et al. (2003) . The open diamond shows the re- 
sult from Tytler et al. (2004). The short-dashed curve 
is the result of Kim et al. (2002) using a sample of 
high-resolution UVES spectra, while the long-dashed 
and dash-dotted curves show the results of Schaye et 
al. (2003) from a combined sample of high-resolution 
UVES and Keck spectra with and without the removal 
of pixels contaminated by associated metal absorption. 
The filled triangles show the results from the LUQAS 
sample for three different redshift ranges (also given in 
Table 2). Damped and sub-damped Lyman-a systems 
have been removed for the estimate from the LUQAS 
sample and errors have been calculated with a jack-knife 
estimate using 40 subsamples. Note that the sample 
used by Schaye et al. (2003) has 14 QSO spectra in 
common with the LUQAS sample. 

The results of PRS and Bernardi et al. (2003), 
which both use low-resolution spectra with comparably 
low S/N, give results which are about 15-20% higher 
than those of studies using high-resolution spectra (see 
also Schaye et al. 2003 for a discussion). As argued 
by Seljak et al. (2003) and Tytler et al. (2004) this 
is most likely due to systematic errors in the contin- 
uum fitting procedure of the low-resolution spectra. We 
agree with this assessment. At redshifts < 2.7 the ex- 
pected continuum fitting errors in high-resolution spec- 
tra is of order 2% (Bob Carswell and Tae-Sun Kim, 
private communication). For spectra at these redshifts, 
the errors are thus almost certainly dominated by the 
somewhat uncertain contribution from metal lines and 
statistical errors. Based on Figure 1 wo will adopt the 
following values for the optical depth and its error in 
our further analysis, Teg{z = 2.125) = 0.17 ± 0.02 and 
Teff (2 = 2.72) = 0.305 ± 0.030. 

3 THE FLUX POWER SPECTRUM OF 
SIMULATED ABSORPTION SPECTRA 

3.1 Numerical code and pEirameters 

We have run a suite of simulations with varying cosmo- 
logical parameters, particle numbers, resolution, box- 
size and thermal histories using a new version of the 
parallel TreeSPH code GADGET (Springel, Yoshida & 
White, 2001). GADGET-2 was used in its TreePM mode 
which speeds up the calculation of long-range gravita- 
tional forces considerably. The simulations were per- 
formed with periodic boundary conditions with an equal 
number of dark matter and gas particles and used the 
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Figure 2. Left panel: ID flux power spectra of different simulations of model B2 at z = 2.75, run with different box-sizes. 
All simulations are with 200^ dark matter particles and 200^ gas particles. Note that simulations have not been scaled to the 
same effective optical depth. The shaded region indicates the range of wavenumbers used to infer the linear dark matter power 
spectrum. Right panel: 3D flux power spectra for different simulations of model B2, for a range of box-sizes and resolutions, as 
indicated in the plot. 



conservative 'entropy-formulation' of SPH proposed by 
Springel & Hernquist (2002). Radiative cooling and 
heating processes were followed using an implementa- 
tion similar to that of Katz et al. (1996) for a primordial 
mix of hydrogen and helium. We have assumed a mean 
UV background produced by quasars as given by Haardt 
& Madau (1996), which leads to reionisation of the Uni- 
verse at « ~ 6, but have also run simulations where we 
artificially increased the heating rates to mimic different 
thermal histories. Most simulations were run with heat- 
ing rates increased by a factor of 3.3 in order to achieve 
temperatures which are close to observed temperatures 
(Schaye et al. 2000; Ricotti et al. 2000; Choudhury, Sri- 
anand & Padmanabhan 2001). 

In order to maximise the speed of the simulations 
we have employed a simplified star-formation criterion 
in the majority of our runs. All gas at densities larger 
than 1000 times the mean density was turned into col- 
lisionless stars. The absorption systems producing the 
Lyman-a forest have small overdensity so this criterion 
has little effect on flux statistics, while speeding up the 
calculation by a factor of ~ 6, because the small dynam- 
ical times that would otherwise arise in the highly over- 
dense gas need not to be followed any more. In a pixel- 
to-pixel comparison with a simulation which adopted 
the full multi-phase star formation model of Springel & 
Hernquist (2003) we explicitly checked for any difi'er- 
ences introduced by this approximation. We found that 
the differences in the flux probability distribution func- 
tion were smaller than 2%, while the differences in the 
flux-power spectrum were smaller than 0.2 %. We have 
also turned off all feedback options of GADGET-2 in our 
simulations. The effect of feedback by galactic winds on 



the statistics of the flux distribution is uncertain but 
is believed to be small (e.g. Theuns et al. 2002, Brus- 
coh et al. 2003, Croft et al. 2003, KoUmeier et al. 2003, 
Desjacques et al. 2004). 

The simulations were all started at z = 99 and we 
have stored 19 redshift outputs for each run, mainly in 
the redshift range 1.5 < z < 3.5. The initial gas temper- 
ature was T = 227 K, and 40 ± 2 SPH neighbours were 
used to compute physical quantities. The gravitational 
softening was set to 2.5 kpc in comoving units for 
all particles. 

We have run a suite of simulations with cosmo- 
logical parameters close to the values obtained by the 
WMAP team in their analysis of WMAP and other 
data (Spergel et al. 2003), as shown in Table 3. Our 
flducial model is a 'concordance' ACDM model with 
r^om = 0.26, Q.OA = 0.74, ^ob = 0.0463 and Ho = 
72kms-^Mpc-^ (B2 in Table 3). The CDM transfer 
functions of all models have been taken from Eisenstein 
& Hu (1999). 

The simulations were run on COSMOS, a shared- 
memory Altix 3700 with 128 900 MHz Itanium proces- 
sors hosted at the Department of Applied Mathematics 
and Theoretical Physics (Cambridge), and a 64-node 
Beowulf cluster with 64 1.2 GHz Sun processors hosted 
at the Institute of Astronomy (Cambridge). The ma- 
jority of our simulations have been run with 2 x 400^ 
particles in a 60 Mpc box and they took about 280 
hrs on 32 processors to reach z = 2. 
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3.2 The effect of resolution and box size on 
the flux power spectrum 

As discussed in K04 and C02 and in Section 2.3, tlie 
range of wavenumbers 0.003 < fc/(s/km) < 0.03 is least 
affected by systematic uncertainties. We will thus limit 
our attempts to infer the linear dark matter power spec- 
trum to this range, which is indicated by the shaded re- 
gion in Figure 2 and following figures. In the left panel 
of Figure 2, we plot the ID flux power spectrum for 
a suite of exploratory runs with fixed particle num- 
ber 2 X 200'^ for five different box sizes/resolutions. All 
simulations are for our fiducial cosmology B2 and have 
been run with the same phases. There is generally good 
agreement in the overlap regions of the flux power spec- 
trum. At large wavenumbers, the flux power spectrum 
for simulations with box size > 60/i^^Mpc, where the 
gas particle mass is larger than 5 x 10^ /i~^Mo, is how- 
ever clearly affected by insufficient resolution. 

The right panel of Figure 2 shows the effect of in- 
creasing the mass resolution by a factor of eight for the 
60/i~^Mpc and 30/i~^Mpc boxes, corresponding to an 
increase of the particle number to 2 x 400^ . Plotted is the 
3D power spectrum which we will use to infer the linear 
dark matter power spectrum using the 'effective bias' 
method developed by C02. There is good agreement be- 
tween the fiux power spectra for k > 0.008 s/km, except 
at small scales for the 60 Mpc box with only 2 x 200^ 
particles. A box size of 60 Mpc with 2 x 400^ parti- 
cles appears to be a suitable compromise for estimates 
of the DM power spectrum. It has converged on small 
scales and will only be moderately affected by cosmic 
variance on the largest scales (for more discussion see 
Section 4.2 on errors). For this choice, the mass per dark 
matter particle is 2 x 10* h~^MQ and the mass per gas 
particle is 4.3 x 10^ /i~^Mq. Note that this is a signifi- 
cant improvement compared to the simulations of C02 
who used 10 dark matter simulations with a box size of 
27.77 ft"^ Mpc and 160* particles. 

3.3 The effect of temperature on the flux 
power spectrum and the rescaling of the 
temperature-density relation 

It is well established by analytical arguments and nu- 
merical simulations that the gas responsible for the 
Lyman-a forest is in photoionisation equilibrium and 
exhibits a rather tight relation between density and 
temperature (e.g. Hui & Gnedin 1997). This relation 
is a quasi-equilibrium established by the balance be- 
tween photo-heating and adiabatic cooling. For moder- 
ate over-density the relation can be approximated by a 
power-law of the form, 

The slope 7 and normalisation Tq of this relation de- 
pend on the reionisation history of the Universe. The 
effect on the flux power spectrum is twofold. Changing 
the temperature changes the width of the absorption 
features and it thus changes the mean flux decrement 
for a given distribution of neutral hydrogen. 



To see how the slope affects the flux power spec- 
trum, it is helpful to consider the fluctuating Gunn- 
Peterson approximation. The relation between the op- 
tical depth and matter density in redshift space can be 
written as 

where f3 = 2.7 — 0.7 7 depends on the temperature- 
density relation of the gas due to the temperature de- 
pendence of the recombination coefficient (see Weinberg 
et al. 1999 for a review). The factor A depends on red- 
shift, baryon density, temperature at the mean density, 
Hubble constant and photoionisation rate. Increasing 
the slope of the temperature-density relation therefore 
flattens the slope of the relation between optical depth 
and matter density. At fixed mean optical depth, this 
results in a decrease of the amplitude of the flux power 
spectrum. 

With a running time of two weeks for each of our 
simulations we could not afford to run an extensive pa- 
rameter study of simulations with different thermal his- 
tories and temperature-density relations. We have thus 
used the rescaling method developed by Theuns et al. 
(1998) to impose a variety of temperature-density re- 
lations. This will not account for the effects that the 
corresponding change in the gas pressure would have 
had on the gas distribution. It nevertheless mimics the 
major effects of different density temperature relations 
on the flux power spectrum reasonably well. To verify 
this explicitly, we have compared two simulations run 
with a factor ^ 3.3 different photo-heating rates. The 
temperatures were different by a factor ~ 2.5. The best 
fitting parameters for the temperature-density relation 
of the 'hot' simulation are To = 10*-^^ K and 7 = 1.6. 

In the left panel of Figure 4, we compare the 
fiux power spectra of the 'hot' and 'cold' simulations 
with flux power spectra for which the temperatures of 
the hot/cold simulations have been rescaled such that 
they have the same temperature-density relation as the 
cold/hot simulations. The differences between the ID 
flux power spectra of rescaled and simulated models 
are smaller than 10% at small k, but start to diverge 
strongly at fe > 0.02 s/km. The differences mostly de- 
pend only weakly on k and the differences between the 
3D spectra will generally be smaller. The middle panel 
shows the effect of rescaling to temperature-density re- 
lations with different 7. We will explore the effect of 
changing the temperature-density relation on the in- 
ferred dark matter power spectrum in more detail in 
Section 4. 

3.4 The amplitude of the matter power 

spectrum and rescaling of the redshift 

Our grid of simulations is somewhat sparse in the fluc- 
tuation amplitude of the matter power spectrum. How- 
ever, at the redshifts considered here, redshift and fluc- 
tuation amplitude are largely degenerate and a suitably 
rescaled simulation output from a different redshift can 
mimic a simulation with different fluctuation amplitude. 
In the right panel of Figure 4 we test how well this 
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Figure 3. Effect of rescaling simulations (60 Mpc box 2 X 400'' particles), to reproduce different temperature-density 

relations and different values of erg. Left: Difference between the ID flux power spectrum of the 'hot' simulation and the rescaled 
'cold' simulation (solid curve) and vice versa (dashed curve). Middle: Difference between the ID flux power spectrum of the 
simulation of the B2 model and the flux power spectra for three rescaled models with different values of the exponent 7 of 
the temperature-density relation. The difference between the ID flux power spectra of two simulations of model B2 with two 
different resolutions is also shown (30 Mpc box, 2 X 400^ particles). Right: The solid curve shows the difference between the 
ID flux power spectra of a simulation of model B2 (erg = 0.85, z = 2.72) and an output of a simulation of model B3 (erg = 1) 
at higher redshift which has the same rms density fluctuation (rescaled to the effective optical depth at z = 2.72). The dashed 
curve shows the difference between the ID flux power spectra of a simulation of model C3 (erg = 1.0, z = 2.125) and an output 
of a simulation of model C2 (erg = 0.85) at lower redshift which has the same rms density fluctuation (rescaled to the effective 
optical depth at z = 2.125). The dotted curve shows the difference between the ID flux power spectra of a simulation of model 
B2 for two different values of r^g. The shaded regions indicate the range of wavenumbers used in our analysis. 



works by comparing the flux power spectrum of simula- 
tions with different ag at redshifts where the rms density 
fluctuation amplitude of the simulation is the same. The 
differences at large scales (fc < 0.03 s/km) are typically 
of the order of 5%. 



4 ESTIMATING THE MATTER POWER 
SPECTRUM 

4.1 The Method 

The method we use to infer the linear dark matter power 
spectrum has been proposed by C02. It uses numeri- 
cal simulations to calibrate the relation between flux 
power spectrum and matter power spectrum. It then 
assumes that the flux power spectrum Ppik) at a given 
wavenumber k depends linearly on the linear real space 
matter power spectrum P{k) at the same wavenumber 
and that both can be related by a simple bias function 
fe(fc), 

Ppik) ^ b\k) Pn,^t{k) . (4) 

The flux power spectrum obtained from the hydro- 
simulations is used to determine b{k). In reality, some 
mode coupling is expected (Gnedin & Hamilton 2002) 



and the relation between flux and matter power spec- 
trum will in general not be linear. However, if the true 
matter power spectrum is close to the matter power 
spectrum used to determine fe(fc) the approximation has 
shown to give accurate results. We refer to C02, Gnedin 
& Hamilton (2002) and Zaldarriaga, Scoccimarro & Hui 
(2003) for more details and for an extensive discussion 
of the possible limitation of this method. 

With our grid of models of varying slope and am- 
plitude of the DM power spectrum, we always have a 
simulation which comes close to reproducing the ob- 
served flux power spectrum to within 10% in the 
range 0.003 < k (s/km) < 0.03. We then use the b{k) 
determined from this simulation to infer the 'observed' 
linear matter power spectrum by dividing the observed 
flux power spectrum by b^{k). The 'corrections' of the 
linear power spectrum with respect to the one actually 
simulated are smaller than 10%. 



4.2 Statistical errors 

As discussed in Section 2.4, the systematic uncertain- 
ties on the observed flux power spectrum are small for 
0.003 s/km < k < 0.03 s/km, the range of wavenumber 
we have used for our analysis. The errors of the ob- 
served flux power spectrum in this range of wavenum- 
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Figure 4. ba^/b for three different simulations with different values of as- Left: ba^/b as a function of T^ff (no scaling of the 
temperature-density relation has been adopted here) . Wc also show the scaling found by Croft et al. (2002) [C02] and Gnedin & 
Hamilton (2002) [GH02] as thin solid and dashed curves, respectively. Middle: ba^/b as a function of Tq, the temperature at the 
mean density, for a fixed Toff = 0.305 and 7 = 1.2. Right: bad/b as a function of 7, the power-law index of the temperature-density 
relation, for a fixed Tea = 0.305 and To = lO'* !^ K. 



bers should be dominated by statistical errors. In the 

following, we will assume the statistical errors of the in- 
ferred linear matter power spectrum to be the same as 
the statistical errors of the observed flux power spec- 
trum. Cosmic variance in the flux power spectrum of 
our simulated spectra is a further possible source of 
statistical error in our modelling. However, our anal- 
ysis is performed using simulations with a box size of 
60h~^ Mpc which corresponds to 6082.3 km/s in veloc- 
ity space. Our simulation thus probes a volume which 
is ~ 27 times larger than (27r/fcmin)^, where femin = 
0.003 s/km is the smallest wavcnumbcr wc use for our 
analysis. The cosmic variance error due to the finite size 
of our simulation box should thus only moderately af- 
fect the smallest wavenumber we use. Its contribution 
to the total error of the rms fluctuation amplitude will 
be negligible compared to the other systematic errors 
which we describe in the next section. 



cal depth increases the flux power spectrum by a factor 

which is nearly constant for k < 0.03 s/km. 

To quantify the efi'ect of the efi'ective optical depth 
on the inferred linear matter power spectrum, we con- 
sider the ratio baa/b, where feed is the average value of 
the bias function b{k) for a fiducial model in the range 
considered, and b is the average value for a model with 
difi'crent parameters (Tett, To and 7). In Figure 4 (left 
panel), we show the effect of changing the assumed ef- 
fective optical depth on the 3D flux power spectrum at 
small k, for three different simulations with three differ- 
ent values of (78 at z = 2.75. The dependence on erg is 
weak. The dependence of bm/b on Teff is well fitted by 

b V0.305y ^ ' 

This dependence is intermediate between that found by 
C02 and that by Gnedin & Hamihon (2002). 



4.3 Systematic errors 

4-3.1 The effective optical depth 

The inferred amplitude of the matter power spectrum 
depends strongly on the assumed effective optical depth 
of the absorption spectrum (e.g. Croft et. al. 1998; C02; 
Gnedin & Hamihon 2002; Scljak, McDonald & Makarov 
2003). In Figure 4 (right panel), we show the effect of 
changing the mean optical depth on the ID flux power 
spectrum. The dotted curve shows the difference of the 
ID flux power spectrum if Teff is changed from 0.305 to 
0.349 for the B2 simulation. Increasing the mean opti- 



4-3.2 The temperature-density relation 

The solid, dotted and dashed curves in Figure 3b show 
the effect of changing 7 on the ID flux power spectrum. 
Increasing 7 leads to a decrease of the flux power, again 
by a factor which is nearly constant for k < 0.03 s/km. 
In the middle and right panels of Figure 4, we show the 
dependence of bfid/b on 7 and To (averaged over k for 
0.003 < fe/(s/km) < 0.03) where baa is again the ratio 
of flux to matter power spectrum for our fiducial model. 
We fixed Teff — 0.305 and rescaled to 7 = 1.2 (middle 
panel) and To = 10'''^^ K (right panel), respectively. 
Rescaled 'hot' simulations were used in all cases. 

The dependences of bfid/b on To and 7 are fitted by 
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b 

and 

fcfid 
b 



i + ro/(io4i5 K) 



1+7/1-3 



(6) 



(7) 



respectively. C02 found a somewhat weaker dependence 
on the value of To, and no dependence on 7. 



4-3.3 Other systematic errors 

Other possible systematic errors include fluctuations of 
the UV background and potentially modifications of the 
forest by galactic winds. Variations of the flux of hydro- 
gen ionising photons are expected to be largo in the 
intermediate aftermath of reionisation. However, at the 
redshifts we are interested in here, they are expected 
to be small due to the large mean free path of ionising 
photons. Meiksin & White (2003) estimate that their ef- 
fect on the flux power spectrum should be smaller than 
0.5% at z = 2.75. 

The effect of galactic winds on the flux power spec- 
trum will depend strongly on the volume filling factor 
of galactic outflows and on when the winds occur (The- 
uns et al. 2002, Bruscoli ct al. 2002). In the range of 
wavenumbers which we used for our analysis. Croft et 
al. (2003) and Desjaques et al. (2004) found a negligi- 
ble effect on the flux power spectrum for galactic winds 
consistent with the correlation between flux decrement 
and galaxies detected by Adelberger et al. (2003) at 
z ^ 3, while Weinberg et al. (2003) found a decrease of 
the flux power at large scales when considering strong 
winds models. Galactic winds appear to affect the flux 
power spectrum mainly by their effect on the strong ab- 
sorption systems which contribute significantly to the 
flux power spectrum (Viel et al. 2004a). 



5 THE INFERRED MATTER POWER 

SPECTRUM 

5.1 The linear power spectrum from LUQAS 
and Croft et al. 

In this Section, we present the main result of the paper, 
the linear dark matter power spectrum we infer for the 
LUQAS and C02 samples. As discussed in Section 4.1, 
we first have to find the simulation output for which 
the 3D fiux power spectrum fits the observed fiux power 
spectrum best. To this end we have used minimisa- 
tion in the range 0.003 < k (s/km) < 0.03. For these 
wavenumbers the covariance matrix is reasonably close 
to diagonal, allowing us to neglect correlations between 
the data points. 

In Figure 5, we show the flux power spectrum of 
the best fitting simulation and compare to the observed 
flux power spectrum for the LUQAS sample (rescaled 
simulation C3 with as — 1.04) and for the C02 sample 
(rescaled simulation C3 with as = 1.035). As discussed 
in Section 2.4, the assumed effective optical depths were 
Teff = 0.17 for the LUQAS sample and Tea = 0.305 



for the C02 sample. We then determined the bias func- 
tion b{k) from the best-fitting models and use Eqn. (4) 
to infer the linear dark matter power spectrum. Using 
the rescaling method of simulations at different red- 
shifts (see Section 3.4), wo obtained a sufficiently fine 
grid in fluctuation amplitude of the matter power spec- 
trum. There was thus no need to interpolate the flux 
power spectrum between simulations of different fluc- 
tuation amplitude. Table 4 gives the inferred 3D linear 
dark matter power spectrum obtained with b(k) from 
the best-fitting simulation in Figure 5 for the LUQAS 
and Croft et al. sample, respectively. The solid triangles 
in the left and right panels of Figure 6 show the cor- 
responding A'iaat = -Pmat(A;) /c^/(27r^). Tlic crrors given 
are the statistical crrors of the flux power spectrum. We 
will discuss systematic errors in the next section. Our es- 
timate of the DM power spectrum from the C02 sample 
is ~ 45% higher than the original result of C02, which is 
shown as the open triangles. The original result of C02 
was obtained with a signiflcantly larger effective opti- 
cal depth of Teff = 0.349 instead of Tes = 0.305. The 
open diamonds show our re-analysis of the C02 data 
with the same effective optical depth of TcH — 0.349 as 
C02 have used. There is good agreement which is re- 
markable considering the fact that we have used hydro- 
dynamical simulations rather than DM simulations and 
that the cosmological model is significantly different 
(f2om = 0.4 VS. f^oni = 0.26). Note, however, that this 
agreement is somewhat fortuitous given that we have 
derived a significantly different scaling of b{k) with Teff 
firom our simulations than C02 do. For our preferred 
value of Teff = 0.305, C02 would have obtained a 25% 
larger amplitude of the power spectrum than we do. The 
solid curves represent a model with n= l,as = 1. The 
dashed curve is the "best fitting" model with a running 
spectral index found by the WMAP team for a com- 
bination of CMB, galaxy survey and Lyman-a forest 
data (Spergel et al. 2003). This model falls significantly 
below the DM power spectrum which we have inferred 
from the LUQAS and the C02 sample. 

5.2 The error budget 

In Table 5, we give a summary of the different sources 
of errors for the inferred rms fiuctuation amplitude of 
the matter density. The statistical error given in Ta- 
ble 5 is that due to the statistical errors of the observed 
flux power spectrum. The systematic uncertainties due 
to Toff, 7 and To are taken from Figure 4. The esti- 
mate of the systematic uncertainty due to the method 
is based on Figure 3. Note that due the weak depen- 
dence of the inferred amplitude on temerature which we 
find the wide adopted possible temperature range nev- 
ertheless leads to a small error due to the uncertainty 
in To. 

The uncertainty due to numerical simulations is dif- 
ficult to estimate, but we believe that we have demon- 
strated that we have sufficient resolution to produce 
convergent results. The dot-daslicd curve in the mid- 
dle panel of Figure 3 shows the effect on the fiux power 
spectrum of a further increase of the mass resolution by 
a factor of eight (this curve can be compared with the 
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0.01 0.01 

k(s/km) k(s/km) 

Figure 5. Left: Flux power spectrum at z = 2.125. The filled triangles are for the LUQAS subsample with a mean redshift 
oi z = 2.125 (see Table 1). The continuous curve is the power spectrum of the simulation that fits the data best. Right: Flux 
power spectrum at 2 = 2.72. The filled triangles are for the fiducial sample of Croft et al. (2002). The continuous curve is the 
power spectrum of the simulation that fits the data best (with t^s = 0.305). Shaded regions indicate the range of wavenumbers 
used in our analysis. 




0.01 0.01 

k(s/km) k(s/km) 

Figure 6. Left: Linear power spectrum at z = 2.125. The linear dark matter power spectrum inferred from the LUQAS fiducial 
sample with the bias parameter derived from a simulation of model C3 which fits the data best. The statistical errors are the 
same as those of the flux power spectrum in Fig. 5. The curve is a theoretical linear power spectrum with parameters as given 
on the plot. The error bar at the bottom of the plot is our estimate of the systematic uncertainty as discussed in Section 5.2. 
Shaded regions indicate the range of wavenumbers used in our analysis. Right: Linear power spectrum at z = 2.72. Empty 
triangles show the result obtained by C02. Empty diamonds are the result of our analysis for the same optical depth as used by 
Croft et al. (rgff = 0.349). Solid triangles are the result of our analysis for our preferred effective optical depth of r^ff = 0.305. 
The solid and dashed curves are theoretical linear power spectra with parameters as indicated in the plot. The dashed curve is 
the "best fitting" model with a running spectral index found by the WMAP team (Spergel et al. 2003). 
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Table 4. 3D linear dark matter power spectrum inferred 
from the LUQAS and Croft et al. (2002) sample. The Croft 
et al. sample was analysed with Teft = 0.305. The values are 
for 7 = 1.6 and Tq = lO^-^^K. 



k (s/km) 



0.00199 
0.00259 
0.00336 
0.00436 
0.00567 
0.00736 
0.00956 
0.01242 
0.01614 
0.02097 
0.02724 
0.03538 
0.04597 



t(A;)(km/s)3 



z=2.125 



z=2.72 



(7.74 ± 120) 
(5.64 ±75.1) 
(3.56 ±46.7) 
(1.68 ±0.62) 
(1.23 ±0.61) 
(4.31 ± 1.71) 
(1.46 ±0.61) 
(9.64 ±2.48) 
(6.77 ±0.59) 
(3.71 ± 0.40) 
(1.78 ± 0.21) 
(9.02 ± 0.98) 
(5.25 ± 0.40) 



X lO'^ 
X lO'' 
X lO'^ 
X 10* 
X 10* 
X 10^ 
X 10^ 
X 10"^ 
X 10^ 
X 10^' 
X 10^ 
X 10^ 
X 10^ 



(9.14 ± 11 
(7.67 ±4, 
(4.61 ± 1, 
(7.52 ± 5, 
(2.77 ± 1.. 
(4.07 ± 0. 
(1.62 ±0, 
(7.09 ±1 
(4.24 ± 0, 
(2.22 ±0, 
(1.05 ±0, 
(7.29 ±0, 
(4.44 ± 0, 



67) X 10* 

47) X 10* 
98) X 10* 
06) X 10^ 
84) X lO'' 
76) X lO'' 
29) X 10''' 
03) X 10" 
75) X 10** 
31) X 10** 
91) X 10** 

48) X 10^ 
21) X 10^ 



Z=2 72 ^ C02x,j^.305 

- - A C02 x,jf=0.349 

§ C02 original 

♦ Seljak et al. (2003) 




-2.90 -2.80 -2. 70 -2. 60 -2.50 -2.40 -2.30 -2.20 

V 



Table 5. Error budget for the determination of the rms fluc- 
tuation amplitude of the matter density field. 



statistical error 


4% 


systematic errors 




Teff(z = 2.125) = 0.17 ± 0.02 


8% 


reff(2 = 2.72) = 0.305 ± 0.030 


7% 


7 = 1.3 ±0.3 


4% 


To = 15000 K ± 10000 K 


3% 


method 


5% 


numerical simulations 


8% (?) 


further systematic errors 


5% (?) 



similar plots in Figure 2). For the relevant wavenunibers 
the difference is less than 2%. However, the discrepancy 
of the scaling of b{k) with Tes compared to C02 is worry- 
ing. We have thus assigned half of the difference between 
their and our result for t^s — 0.305 as systematic uncer- 
tainty due to numerical simulations. This is admittedly 
arbitrary and merits further investigation. 

Further unknown systematic uncertainties are ob- 
viously impossible to quantify. We have nominally as- 
signed 5% to take into account a possible effect due to 
galactic winds. This could be larger especially at low 
redshift for the LUQAS sample, where less is known 
observationally about the effect of galactic winds. The 
sum of the systematic errors (added in quadrature) for 
the fluctuation amplitude is 14.8% at a = 2.125 and 
14.3% at z = 2.72, and is shown as the error bars at the 
bottom of Figure 6. 

5.3 Comparison with previous estimates 

In Figure 7, we compare the constraints on the ampli- 
tude and slope of the power spectrum with previous 
estimates. We fit the linear dark matter power spec- 
trum with a power-law function P{k) = Pp{k/kpY with 
kp = 0.03 s/km, as in C02. We used a diagonal likelihood 



Figure 7. Constraints on the slope and amplitude of the 
linear dark matter power spectrum for the fiducial sample of 
C02 {z = 2.72). The solid contours show the 68.3% and 95.4% 

confidence levels for our analysis with T^g = 0.305, while the 
dashed contours are for r^g = 0.349. The empty diamond 
shows the original result obtained by C02 with la error bars. 
The filled diamond shows the best fitting values obtained by 
Seljak, McDonald & Makarov (2003) (for t^s = 0.30). 



for this estimate and convolved the likelihood in the Pp 
direction with a Gaussian function to take the system- 
atic errors into account (eqn. 15 of C02). The contour 
levels show the 68% and 95% confidence levels in the 
amplitude-slope plane for Tf.« = 0.305 (continuous line) 
and Tefi = 0.349 (dashed line), respectively. 

The empty diamond with error bars is the determi- 
nation by C02, which is in good agreement with our de- 
termination for the same Tes- As expected, the inferred 
amplitude incrcEises significantly if the assumed t^s is 
reduced. The result of Seljak, McDonald & Makarov 
(2003) for ToEF = 0.3 (their Fig. 1) is also shown as the 
filled diamond and agrees well with our determination. 
Note, however, that the Seljak et al. result has been 
obtained differently. Seljak et al. fitted a large grid of 
flux power spectra obtained from HPM (Hydro-Particle- 
Mesh, sec Gnedin & Hui 1998; Meiksin & White 2001) 
simulations with six free parameters to the observed flux 
power spectrum of the C02 sample. 



5.4 Combining CMB and Lyman-a forest data 
to constrain n and erg 

The measurement of the amplitude of the matter power 
spectrum on scales of a few Mpc with the Lyman-a for- 
est is a powerful tool to constrain the spectral index 
of primordial density fluctuations when combined with 
measurements on large scales from CMB fluctuations 
(see Phillips et al. 2001 for a detailed discussion). The 
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^eff '^eff 

Figure 8. Left: The values of and n of the best fitting COBE normaUzed Hnear power spectra as a function of the assumed 
effective optical depth, for the linear power spectrum inferred from LUQAS at z = 2.125. The solid curve is obtained with 
numerical simulations with 7 1.6 while the dashed curve is for simulations which were rescaled to have a density temperature 
relation with 7 = 1.0. Shaded regions indicate our preferred values for the effective optical depth as obtained from high-resolution 
spectra. The arrow indicates the estimate obtained by Press, Rybicki & Schneider (1993) from low-resolution spectra. Right: 
The same for our re-analysis of the Croft et al. (2002) data at z = 2.72. The dotted line in the right panel is the lower redshift 
result scaled by a factor 0.305/0.17 in the x-axis. 



linear matter power spectrum at 2; = can be written 

as, 

PnMk)^Ak"T^{k), (8) 

where T{k) is the matter transfer function which de- 
pends on cosmological parameters in the usual way. 

To get a first idea how the spectral index inferred 
by such a combined analysis depends on t^s, we here as- 
sume a COBE normalized power spectrum and no con- 
tribution by tensor fluctuations. We then find the best 
fitting spectral index n (having fixed the other cosmo- 
logical parameters to our fiducial values) . We use CMB- 
FAST to calculate the theoretical linear power spectra 
(Seljak & Zaldarriaga 1996). Note that normalizing to 
the WMAP data would give similar results. The thick 
solid curve in Figure 8 shows the result for a range of 
Teff using numerical simulations which have 7 ~ 1.6 and 
To fa 10*'^^ K. The left and right panel are again for 
the LUQAS and C02 sample, respectively. The shaded 
regions indicate the range of preferred by high- 
resolution absorption spectra as discussed in Section 2.4. 
For the errors we have added the statistical errors and 
the systematic errors in Table 5 in quadrature. Note 
that in Figure 8 we have omitted the error due to TbB 
and 7 as the dependence on these parameters is shown 
explicitely. 

Once the spectral shape and the cosmological pa- 
rameters are fixed, the measurement of the fluctuation 
amplitude on Lyman-a forest scales can also be ex- 
pressed in terms of as , the rms fluctuation amplitude of 
the density for a 8 Mpc sphere. The corresponding 



values are shown on the right axis of Figure 8. Note that 
the scale corresponding to 8 Mpc is ~ 0.008 s/km 
ai z — 2.72 (assuming again the cosmological param- 
eters of Section 3). The flux power spectrum is thus a 
direct probe of ag for the range of fc-values used in our 
analysis. 

The inferred values of as and n will also depend 
on the assumed cosmological parameters (see Phillips 
et al. 2001 for a detailed discussion). Most important 
is the dependence on Qom and h. For a flat cosmolog- 
ical model, the spectral index scales approximately as 
n <x (r^om/i^)""'^^ (Phillips et al. 2001). In Figure 9, we 
show the depencence of the inferred value of as on Qom 
for a flat universe. This dependence is essentially negli- 
gible. The dependence on h is also weak. For the error 
estimate we have added all values in Table 5 in quadra- 
ture. Table 6 lists our final estimate of as and n for 
the LUQAS and C02 sample, as well as their weighted 
mean. 

5.5 Gravitational growth 

Our subsample drawn from the LUQAS sample was cho- 
sen to maximise the contrast in redshift with respect to 
the C02 sample, and to further investigate the redshift 
evolution of the flux power spectrum. There are two ef- 
fects which are responsible for the evolution of the flux 
power spectrum, the decrease of TbB with decreasing red- 
shift, and the increase of the fluctuation amplitude of 
the matter power spectrum due to gravitational growth. 
As can be seen by comparing the right and left panels 
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Table 6. Estimated values of as and n and their errors (7 = 1.3 and Teff = 0.17 and t^s = 0.305 a,t z = 2.125 and at 2 = 2.72. 
Note that n scales as n oc (nom/i^/0.135)-°-35 (a). 



LUQAS 


C02 


combined 


0-8 0.95 ± 0.04(stat.) ± 0.13{syst.) 
n 1.02 ± 0.02(stat.) ± 0.08(syst.) 


0.92 ± 0.04{stat.) ± 0.13(syst.) 
0.99 ± 0.02(stat.) ± 0.08(syst.) 


0.93 ± 0.03(stat.) ± 0.09(syst.) 
1.01 ± 0.02(stat.) ± 0.06(syst.) 



(a) Note that there is an additional uncertainty in the COBE normalization (Bunn & White 1997). 



1.4 



1.2 



^1.0 



0.8 



CROFT z = 2.72 
LUQAS = 2.125 



0.6 



0.20 



0.25 



0.30 



0.35 
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Figure 9. The value of as of the best fitting COBE nor- 
malized linear power spectra as a function of the assumed 
fiom for the linear power spectrum, inferred from LUQAS at 
2 = 2.125, and from the C02 data at 2 = 2.72. The assumed 
effective optical depths are t^s = 0.17 and t^s = 0.305 
and the exponent of the temperature-density relations was 
7 = 1.3. 



of Figure 5, the net effect is a decrease of the flux power 
with decreasing rcdshift which is, however, significantly 
smaller than that expected from the decrease of t^b- 

In order to assess if the flux power spectra of the 
LUQAS and C02 samples are consistent with the ex- 
pected gravitational growth of the matter power spec- 
trum between the two redshifts, we can compare the 
inferred values of cts, which should then be the same. 
To facilitate such a comparison, the dashed curve in the 
right panel of Figure 8 is the inferred as of the LUQAS 
sample (shown in the left panel) with all effective opti- 
cal depths scaled by the same factor of 0.305/0.17. The 
inferred values of as agree to within the errors, and the 
evolution of the flux power spectrum between z = 2.7 
and z = 2.1 is fully consistent with being due to the ex- 
pected gravitational growth and the observed evolution 
of Teff. The same was found by C02 when they compared 
their fiducial sample to their low-redshift subsample. 



6 DISCUSSION AND CONCLUSIONS 

We have used the observed fiux power spectrum of two 
large samples of high- resolution spectra, a sample drawn 
from LUQAS at a median redshift of 2; = 2.125, and the 
sample compiled by Croft et al. (2002) at = 2.72, 
together with a suite of high-resolution numerical sim- 
ulations, to infer the dark matter power spectrum on 
scales 0.003 < fc(s/km)< 0.03. 

We have obtained the following results: 

(i) With the same assumptions for effective optical 
depth, density-temperature relation, and cosmology, our 
inferred linear matter power spectrum agrees very well 
with that inferred by Croft et al. (2002). 

(ii) We confirm previous results that the inferred rms 
amplitude of density fluctuations depends strongly on 
the assumed Teff . It increases by 20% if we assume an 
optical depth of TcS ~ 0.305 a value suggested by studies 
of high-resolution absorption spectra. We find, however, 
a dependence on Teff which is weaker than that of C02 
and stronger than that of GII02. 

(iii) For values Teff suggested by high-resolution ab- 
sorption spectra the linear power spectrum of the best 
fitting running spectral index model of Spergel et al. 
(2003) falls significantly below the linear power spec- 
trum inferred from both the LUQAS and the C02 sam- 
ple. 

(iv) The decrease of the amplitude of the fiux power 
spectrum between z = 2.7 and z = 2.1 is consistent 
with that expected due to the decrease of Tea and the 
increase of the amplitude of matter power spectrum due 
to gravitational growth. 

(v) Our estimate of the systematic uncertainty of the 
rms fluctuation amplitude of the density (~ 14.5%) is 
a factor 3.5 larger than our estimate of the statistical 
error (~ 4%). The systematic uncertainty is dominated 
by the uncertainty in the mean effective optical depth 
and — somewhat surprisingly — by the uncertainties 
between the numerical simulations of different authors. 
Reducing the overall errors will thus mainly rely on a 
better understanding of a range of systematic uncertain- 
ties. 

(vi) By combining the CMB constraint (assuming 
that there is no contribution from tensor fluctua- 
tions) on the amplitude of the DM power spectrum 
on large scale with the high-resolution Lyman-a for- 
est data we obtain n = 1.01 (nom/i^/0.135)"°'^^ ± 
0.02(statistical) ± 0.06(systcmatic) for the spectral in- 
dex. The corresponding rms fluctuation amplitude is, 
as = 0.93 ± 0.03(statistical) ± 0.09(systematic). 
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